Optomechanical back-action evading measurement without parametric instabihty 
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We review a scheme for performing a back-action evading measurement of one mechanical quadra- 
ture in an optomechanical setup. The experimental application of this scheme has been limited by 
parametric instabilities caused in general by a slight dependence of the mechanical frequency on 
the electromagnetic energy in the cavity. We find that a simple modification to the optical drive 
can effectively eliminate the parametric instability even at high intracavity power, allowing realistic 
devices to achieve sub-zero-point uncertainties in the measured quadrature. 



I. INTRODUCTION 

In combination with back-action evading techniques 
described below, squeezed states offer the potential for 
important applications in optomechanical precision force 
sensing [l| , in particular achieving sensitivity with resolu- 
tion below the zero-point fluctuations Xzp of the mechan- 
ical component of the sensor. Unfortunately, due to the 
typical thermal environment of a mechanical system, it is 
difficult to produce direct squeezing below the zero-point 
level in a mechanical oscillator. It is similarly challenging 
to ascertain its position to that level of accuracy in a sin- 
gle strong measurement. One possible solution is to con- 
tinuoush' observe the position in, e.g., an optomechanical 
setup [^, using light transmitted through a Fabry-Pcrot 
cavity to probe the motion of an oscillating end-mirror or 
using an equivalent microwave circuit with mechanically 
modulated frequency. Such a measurement process grad- 
ually reduces Ax. However, this leads to a problem as 
the oscillator rotates through phase space. Because mea- 
suring X increases the uncertainty of p and every quarter 
cycle X and p are exchanged, the maximum squeezing is 
limited to roughly that which could be produced in only 
one fourth of the oscillator period, even if technical noise 
is neglected. Hence, the potential squeezing is limited by 
measurement back-action. 

A back-action evading measurement of the position 
of a membrane in a cavity optomechanical system was 
proposed as early as 1980 in Refs. [3, |j], who suggested 
driving the resonator with an input field resonant with 
the cavity frequency Uc, but modulated at the mirror 
frequency f2m- In Fourier space, the field has oscilla- 
tory components at ujc ± f2,„; hence, this scheme is often 
known as two-tone back-action evasion. By modulating 
the light field frequency at Vim the measurement effec- 
tively turns on and off as the system oscillates. This pro- 
tocol thereby measures neither position nor momentum 
individually, but rather, one of the mechanical quadra- 
tures. Thus, while the measurement back-action still ex- 
ists, it feeds only into the unmeasured quadrature and 
evades the measured one, leaving in principle no lower 
limit on the uncertainty one quadrature might reach. 

A detailed quantum mechanical analysis is presented 
in Rcf. [6|, and a generalized version exploiting interef- 



erence between additional tones is proposed in Ref. [7|. 
The technique was originally explored for use in grav- 
ity wave detectors [5|, including an early approximation 
to stroboscopic position measurement [8|, but these ex- 
periments using massive oscillators were not designed to 
reach sensitivities near their zero-point levels. More re- 
cently, experiments have been carried out in the quan- 
tum regime. These have reached sensitivities of Ax^^ [0], 
2.5xzp [l^, and lAx^.^ [il|, but, thus far, no experiment 
has achieved sub-zero-point position sensitivity. 

While two-tone back-action evasion is an elegant so- 
lution, experimental reality intervenes to place a rather 
restrictive limit on such a scheme. Because the enve- 
lope of the driving field oscillates at the mechanical fre- 
quency, the intracavity power oscillates at twice the me- 
chanical frequency. Through indirect effects, this leads 
to the frequency of the mechanical oscillator becoming 
slightly modulated, with the modulation oscillating at 
twice the natural frequency. Such a frequency modula- 
tion produces a parametric instability, which will drive 
the system and greatly reduce the amount of squeezing 
possible. To accomplish the back-action evading mea- 
surement, high optical power and low mechanical dissi- 
pation are both critical, yet these factors both worsen the 
parametric instability. 

The variety of ways in which this instability can arise 
is staggering. In the experiments mentioned above, it 
originated respectively from non-linear terms in the op- 
tomechanical coupling Q , cavity heating causing a ther- 
mal shift in the frequency of the mechanical element [10| , 
and two-level systems in surface oxides acting as non- 
linear dielectrics |ll| . Thus, the common limitation of 
these experiments is a connection between mechanical 
frequency and the cavity energy, 5Q.rn (x. E. In the two- 
tone scheme, this connection seems to lead inevitably to 
parametric instability. 

However, it is in principle possible to work around 
this particular limitation. There are two critical features 
needed to avoid the parametric instability while perform- 
ing the back-action evading measurement: 

1. The power in the cavity must not oscillate at twice 
the mechanical frequency. 

2. The probe light must couple only to a single 
quadrature of motion. 



The first of these requirements prevents the instability, as 
we will discuss in more detail in the next section. The lat- 
ter requirement prevents measurement back-action from 
affecting the measured quadrature. Small deviations 
from this requirement will reduce the ultimate sensitiv- 
ity of the measurement, but do not preclude sensitivities 
below the zero-point level, as shown in Ref. [6|. 

We will verify in this paper that these two criteria may 
be pursued somewhat independently of each other. In the 
next section, we show that the ideal envelope for the in- 
tracavity field is a square wave, which satisfies the first 
criterion above. While experimental constraints prevent 
the realization of a perfect square wave, we show that it 
is sufficient to add a single additional drive tone. The 
rest of the paper is devoted to proving that the second 
criterion will also be satisfied. We review in Section III 
the optomechanical system used for the measurement, in- 
cluding a general optical driving term and outlining the 
main steps in extending the results of Ref. [12l | for such a 
drive. Using reasonable simplifications, in Section IV we 
then derive a general master equation for the conditional 
evolution of the mechanics under continuous observation 
and show in Section V that, for the desired purpose, it 
can be reduced to the master equation of Ref. [y] , thereby 
completing the demonstration of the consistency of the 
two criteria. Finally, Section VI is a summary and out- 
look. 



II. FIELD ENVELOPE 

We model a generic cavity optomechanical system as 
a single-mode optical or electrical resonator with reso- 
nant mode frequency ujc whose radiation pressure drives 
an harmonically confined end-mirror or capacitor plate 
with natural oscillation frequency ^m- The cavity is 
driven on-resonance but modulated with a field envelope 
a[t). In the microwave domain, this cavity field is often 
achieved through the application of multiple tones in or- 
der to avoid the low frequency phase noise of available 
sources. 

We proceed by first asking what specific drive scheme 
can be used to avoid the mechanical parametric instabil- 
ity. To recap, this instability is due to a generic coupling 
arising through a variety of mechanisms between the me- 
chanical frequency and the electromagnetic energy in the 
cavity which oscillates proportional to |a(f)p. If only two 
tones are used to drive the system at Wc ± f^m, the cav- 
ity power then oscillates at 2i7m. Via device-dependent 
mechanisms, this induces a shift in the mechanical fre- 
quency. 



^L(t)=^r, 



5Q,m cos2ri„,t. 



(1) 



In a well-made device, the shift 5Vl,n will be much less 
than Vtjm though the quality factor Q will also be high. 
Hence the mechanical damping will be weak. The onset 
of parametric instability begins when the criterion for the 



fractional frequency shift 
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is satisifed[13{. Similar parametric resonances occur if 
the mechanical frequency oscillates at subharmonics of 
2ilm, but our scheme avoids these. Furthermore, the re- 
quired fractional frequency shift is of order Q~^'" for the 
nth subharmonic, rendering all but the first instability 
irrelevant for our purposes. However, even in carefully 
engineered devices, the limit in Eq. ([5]) is reached at low 
pump power, reducing the maximum possible accuracy of 
the back-action evading measurement to above the zero 
point level|ll|. 

We will show in the following sections that the mea- 
surement of the membrane motion can be performed in 
such a way that it is dominated by the two tones at 
LOc ± ^m- Therefore, we are more or less free to add ad- 
ditional fields to the cavity, as long as they are not near 
those frequencies, without affecting the measurement- 
based squeezing. We can in that way cancel out the 
oscillations of the energy in the cavity, or at least push 
them to higher harmonics, where the instability should 
disappear. 

Consider specifically a field envelope, i.e., the electro- 
magnetic field in a frame rotating at cUc, with a single 
added drive tone at Silm, 

a(t) ex cos n„^t + ^e*(3f2™t-h*) (3) 

The energy E in the cavity is proportional to |ap. 



E{t) ex -+ fi'^ + A2 cos {2fl,nt + $2) + Ai COS {4n,nt + $) , 



with 



(4) 



A2^ \l-+ fiCOS^ + H'^ 



$2 = arctan 



2/1 sin $ 
1 + 2^ COS (f> ' 



(5) 



Thus, by adding the additional 3J7m-detuned tone with 
/i = 1/2 and <I> = tt, the energy oscillating at 2f7,„ is 
completely redirected via interference into the DC and 
4ri„i terms. For small deviations from ideality, i.e., 
fi ~ 1/2 + (5/i, $ = TT + (5$, the amplitude of the en- 
ergy oscillation at 2flm is 



A. 



{s^,y 



(5$)2 



(6) 



which still represents a significant reduction in the 2f2,„ 
Fourier component. The drive needed to produce such 
an envelope can easily be deduced by inverting Eq. p^ 
found in the next section. If desired, one could continue 
to eliminate the higher harmonics of the energy oscilla- 
tion by adding additional phase-locked sources detuned 
from the cavity by ±{2n + l)n„i, for increasing integer 



n. If we extend the series in Eq. ([3]), use tones of equal 
strength at each pair of red- and blue-detuned sidebands 
so that a{t) is real, and repeat the cancellation proce- 
dure of Eq. ^ for these higher harmonics, we find that 
a{t) converges to a square wave envelope - a function 
that oscillates at the same fundamental frequency as the 
mechanics but which, when squared, is simply a con- 
stant. Thus, the energy in the cavity remains constant, 
and there is no resulting oscillatory mechanical frequency 
shift at any harmonic of the mechanical motion. In hind- 
sight, such a solution is self-evident. In terms of the field 
itself rather than the envelope, this corresponds to a tt 
phase shift every half mechanical period. 

However, adding numerous phase-locked tones may be 
technically difficult, and so we examine the effectiveness 
of a single added Srim-detuncd tone. The energy in the 
cavity, and hence the mechanical frequency, now has a 
Fourier component, aside from the DC term, at 4r2,„. In 
this case, the required fractional frequency shift to excite 
a parametric instability is greater than 1 [13| . The power 
necessary to cause a shift of that magnitude far exceeds 
that needed to perform the back-action evading measure- 
ment and would likely have additional deleterious effects 
on the mechanical device. We therefore find that a sin- 
gle additional tone should be sufficient to eliminate the 
instability. 



III. MODEL 

Having now satisfied the first key criterion listed in the 
introduction we turn to the demonstration that, under re- 
alistic conditions, the probe light couples predominantly 
to a single quadrature of motion. Our starting point is 
the generic optomechanical Hamiltonian, see e.g. Ref. Q 

H = h{ujc — Gx) [d^ a — (a^a)(t)] + MlmC*c 

+ ih^{bUi)a-h,,{t)a)), (7) 

that describes the dynamics of a single-mode Fabry-Pcrot 
resonator with an oscillating, harmonically bound end- 
mirror driven by radiation pressure. Here d is the anni- 
hilation operator of the cavity mode and k its decay rate, 
6in is the amplitude of the (classical) driving field, and 
X = a;zp(c-|-c^) is the displacement of the mirror resulting 
from radiation pressure. G is the single photon optome- 
chanical coupling. Finally x^p is the position variance of 
the membrane in its ground state. 

We assume that the drive is resonant with the cavity 
but periodically modulated at the mechanical frequency. 



hn{t) 



E/^^ 



iVk^lt 



(8) 



This differs from the scheme of Refs. [j, Q where N — 1 
only /3-I-1 ^ 0. We further assume that the mechanical 
element is thermally damped at a rate 7 due to its contact 
with a reservoir of temperature T. Thus its equilibrium 



thermal phonon occupancy is given by the Bose-Einstcin 
distribution. 



1 



exp(?ir2„j/fcBT) - 1' 



(9) 



However, we assume that fc^T <C hiUc so that we can 
neglect thermal photons. 

The effects of thermal coupling and the cavity decay 
are encapsulated by the master equation for the compos- 
ite density matrix p of the oscillator-field system. 



dp 
'dt 



7 



[p, H] + -V[d]p +Uin+ mi^P + n-D[c^p) 



where the dissipation super-operator V is defined as 
'D[o\p = 26p6' — o* dp — pd* d. 



(10) 



fll) 



We also work in the good cavity limit, specifically assum- 
ing a separation of scales 



7 <C K <C r^m ^ <^c 



(12) 



We now simplify the master equation by moving into 
the rotating frame for both the optical mode and the me- 
chanical oscillator, and then displacing a by a classical 
mean value a(t) - ultimately, the field envelop of the pre- 
vious section - where the time dependence here is due to 
the modulated nature of the drive. Mathematically, this 
is equivalent to applying several unitary transformations 
to the density matrix in succession 



where 



p' ^WVU p\j'^V'^W\ 

U = exp(irimC^ci), 
V = exp{iu>cd'dt), 
W = cxp[ait)d'' - a* {t)d]. 



(13) 



(14) 



Though W{t) docs not typically commute with W{t'), 
this fact ultimately contributes only an irrelevant net 
global phase to the evolution. Identifying (a'I'd) = |ap 
and noting that in the rotating frame the mean intracav- 
ity field satisfies 



da/dt + Ka/2 = v^^ine 



iujct 



(15) 



we find readily 



ait) 



ae 



00 



a £6 



n„^et 






(16) 



The Hamiltonian governing the evolution of the system 
in the primed frame is 

H'^g (ce-*"-* + cte'^^"*) {a{t)d^ + a* [t)d + 0+0) , 

(17) 



where g = Gx^p. The dissipative terms remain un- 
changed, except for the replacement of p with p'. In 
subsequent calculations, we neglect the d^a term in the 
Hamiltonian (|17p as its effects on the system's evolution 
are significantly smaller than those of the preceding two 
terms, which are multiplied by the classical mean a{t). 



IV. OPTOMECHANICAL MASTER EQUATION 

We now turn to the issue of measurement of the sys - 
tem. Our derivation follows the approach of Ref. [13 |. 
generalized to a local oscillator with a time-dependent 
amplitude or phase. 

By using a homodyne detection scheme, it is possible 
to make a measurement of one quadrature of the cavity 
field. The output field from the cavity is 



Hamiltonian and simplifying it to the form given below 
in Eqs. ([^ and (P5|) . Second, we can adiabatically ehm- 
inate the elements of the density matrix that are off- 
diagonal with respect to the optical field. Finally, this 
will allow us to trace out the light field and derive an 
effective master equation for the mechanical subsystem 
alone. Because similar calculations have been reported 
elsewhere, including in Ref. [l2|, we do not reproduce all 
intermediate details here. 

Our first step is to make the RWA, after which we are 
left with the interaction 



H' 



Cla 



Cio^ 



where 



C*! = aic + a_ic 



(22) 



(23) 



601 it — Oi„ 



(18) 



and 



and the field reaching the detector is -Bio + ^out, where B\o 
is the additional local oscillator. We can fold all the clas- 
sical contributions together into a "net" local oscillator 
strength given in the rotating frame by 



S„ 



(Sio + &i„)e'"=7^/^ - a{t) = S(i)e'^(*' (19) 



where (p is relative phase between the local oscillator and 
the ouput field. Note that if B\o is sufficiently large, the 
additional terms are negligible. For a detector of effi- 
ciency T] the detected photocurrent / (i.e., the measure- 
ment record) is given by 

Idt = {B^ + B{ae-'"^ + a) e''>')) ijndt + By/Tf^dW, (20) 

where W{t) is a Wiener process. That is, ^(t) = dW/dt 
is Gaussian white noise and (dW)'^ = dt. The dynamical 
effects of this measurement on the density matrix can 
be computed and yield the conditional stochastic master 
equation (SME) [l3 



dp. 
dt h"'"'" ' ■ 2 

, 7 



-[p,,H'] + -V[d]p, 



j [{n + l)V[c]pc + ?lP[ct]pc] 



(21) 



+ {ape-"'' + pa)e">' - {de-">' + a^e}'>')p) y^S.{t). 



V. REDUCED MASTER EQUATION AND 
MEASUREMENT 



In order to proceed, we now adiabatically eliminate 
the light field. This can be achieved by working in the 
weak-coupling limit. Specifically, we assume that the op- 
tomechanical interaction frequency, approximately given 
by g\a{t){c)\, is much less than the cavity decay rate k. 
This has multifold advantages: First, we can accurately 
make the rotating wave approximation (RWA), thereby 
removing the explicit time dependence of the interaction 



[Ci,Cl 



"1 



«-i 



(24) 



If the envelope modulating the light field is real, then, for 
all i, a*f, = a-i, and Ci becomes a multiplicative constant 
times a Hermitian mechanical quadrature operator Q. 
Namely, 



Ci = AV2Q, 



(25) 



where we have taken a±i — Ae^^^ , with A real, so that 
the Hamiltonian ultimately reduces to. 



H' = gAV2 (d + dt) Q. 
The next step is to define the small parameter 



gA 



(26) 



(27) 



Because of the dominant effect of dissipation on the dy- 
namics of the intracavity light field it always remain near 
its equilibrium, absent the optomechanical interaction. 
In the displaced frame, this is the ground state. There- 
fore, we can expand the density matrix in the (displaced) 
optical Fock basis as 

P = Poo|0)(0|+6(poi|0)(l|+H.c.) 
+ eVii|l>(l| + e'(P02|0)(2| + H.c.) + ^(e^). (28) 

The adiabatic elimination proceeds by substituting equa- 
tion ((28)) into the stochastic master equation (|2ip . We 
can then obtain a dimensionless version of the SME by 
rcscaling time and the Wiener increment, dr = ndt, dw = 
y/ndW. (Note that dw/dr is still Gaussian white noise 
because {dwY = dr.) Adiabatically eliminating the off- 
diagonal terms and truncating after order e^ yields the 
relations 



P02 = 2ipoiQ, 



Poi 



2iV2pooQ, 



(29) 



and the equations of motion for the diagonal terms arc 
(again, to order e^) 



dpo 



iV2(poiQ~ Qplij +Pii 



e^dr 



poie** + pl,e-'^ - pootr Uie** + pl.e'^^ 



edw. 



dpii = - (i\/2 [Qpai - pliQj + pii dr. 



(30) 



where the trace in Eq. (PU)) is over the mechanical degree 
of freedom and we have omitted the dissipation terms for 
brevity. 

Tracing over the optical degree of freedom is equivalent 
to finding an expression for d{pQo + e^pn), which is done 
easily by substitution of Eq. ^^ into Eqs. ([30]). After 
substituting back in t and W, this yields the stochastic 
master equation for the conditional density matrix of the 
mechanical subsystem, 



dp„ 



dt 



-k[Q,[Q,p,n]] 
V2^{e'^p,nQ 



^Qpm - 2i s'm <f){Q)p,n)£,{t) 



-({n + l)V[c]p,n + nV[c^p,^), 



(31) 



where k ~ Ag'^A^/n. 

The relative phase of ai and a_i defines which quadra- 
ture of motion is measured and without loss of generality 
we can take 6' = so that Q = X\ we use Y to represent 
the orthogonal quadrature. On the other hand, Eq. ([51]) 
emphasizes the importance of the relative phase <f> be- 
tween the local oscillator used for detection and the clas- 
sical driving field. If they are 7r/2 out of phase, then max- 
imum measurement strength is achieved. On the other 
hand, if they are in phase (0 = 0), then the measurement 
acts as a stochastic unitary drive of the system, i.e., a 
random force displacing the mechanics. In this case, the 
light quadrature being measured is d-\-a) , which amounts 
to replacing that term in the Hamiltonian with a fluctu- 
ating, semi-classical value. 

We proceed by taking = it/2, which reproduces the 
master equation of Ref. Q (up to a trivial net minus 
sign due to the difference in sign between Eq. (pS)) and 
the Hamiltonian of that paper) , 



dpn 



dt 



= -k[X, [X, p„,]] - ./2^{prnX + Xprn - 2(X))e(t) 

+ l{{n + l)V[c]p + nV[c^]p). (32) 



The equations for the conditioned mean values, vari- 
ances, and covariancc (C = {XY + YX)/2- {X){Y)) arc 



readily derived under a Gaussian state ensatz: 

^^iX)^-liX)-VKV.^, 

dV: 



X 



dt 

dVy 

"df 



KC^ + 2/c - 7 ( W - n 



^^-KVxC-jC, 



(33) 

(34) 
(35) 
(36) 
(37) 



where we have introduced the scaled measurement 
strength K = Srjk. These variances approach steady 
state values. Of particular interest is the variance of X 



Vx-^|^{2n- 



1+2^) -2^' (^«) 



which approaches for sufficiently large K or small 7; the 
uncertainty relations are maintained by a concommitant 
increase in Vy- These results are now identical to those 
of Ref. [a] • Because our modified system produces all the 
same results from this point, we can omit further cal- 
culations on the use of feedback to promote conditional 
squeezing into unconditioned, or "real" , squeezing. 

Summarizing, then, we have shown that, provided that 
the RWA can be invoked, only the contribution of the first 
sideband, £ = ±1, contributes significantly to the homo- 
dyne detection signal, demonstrating that the probe field 
couples predominantly to a single quadrature of motion. 
With this, both criteria discussed in the introduction 
are satisfied, demonstrated the validity of the proposed 
scheme in eliminating or strongly reducing the impact of 
the parametric instability. 



VI. CONCLUSIONS 

We have thus extended the proposal to perform op- 
tomechanical back-action evading measurements to the 
case of a multi-tone drive scheme. Specifically, by sim- 
ply adding a third tone detuned from the cavity by S^m 
with appropriate amplitude and phase, we push the os- 
cillations in cavity energy to higher harmonics of the me- 
chanical frequency, which in turn do not contribute to the 
parametric instability. We then reproduced, in the ap- 
propriate good-cavity, weak-coupling regime, the desired 
back-action evading measurement of a single mechanical 
quadrature. Because these instabilities appear in such 
diverse experimental settings, we hope this modification 
will prove useful in reaching sub-zero-point position sen- 
sitivities deep in the quantum regime. 

There is still room to extend the results presented 
above. For instance, it may be instructive to include 
the frequency modulating terms explicitly in the Hamil- 
tonian, ii/parametric c)c o' ttc' c. In addition, we can go 
beyond the rotating wave approximation to examine the 



back-action contributed by the rapidly oscillating terms 
in the Hamiltonian. While small, this contribution is 
non-zero, and could eventually impose a lower limit once 
other technical obstacles are overcome. Finally, measure- 
ments of the output field other than homodyne detection 
should be considered to increase the applicability of the 
scheme even further. 
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